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ABSTRACT 

The  aaalyaa  cl  tha  tnoiaat  tmpeum  of  a  limit  laometry  diak  electrode  (embedded  ia  aa  infinite 
insulating  plane)  nihjecied  to  aa  AC  flux  perturbance  is  dircimad  The  real  and  imaginary  components 
of  the  impedance  era  derived.  The  analysis  makes  use  of  the  properties  of  discontinuous  definite  Beeael 
integrals. 


INTRODUCTION 

jj 

We  have  shown  recently  [1,2]  that  it  is  possible  to  derive  the  chronopoteatiomet- 
ric  and  chronoamperometric  response  of  microdisk  electrodes  by  modelling  these 
experiments  with  a  constant  uniform  flux  over  the  microdisk  and  by  deriving  the 
average  concentration  over  the  surface.  The  properties  of  the  discontinuous  integrals 
of  Bessel  functions  are  used  in  these  calculations  to  account  for  the  mixed  boundary 
conditions  in  the  plane  of  the  electrode  (prescribed  flux  over  the  electrode  surface, 
zero  flux  over  the  surrounding  insulator  surface).  It  is  evident  that  the  use  of 
disontinuous  integrals  allows  the  discussion  of  a  wide  range  of  electrochemical 
techniques  and  in  this  paper  we  discuss  the  derivation  of  the  AC  impedance  of  a 
microdisk  electrode  for  a  simple  redox  reaction.  f.  t  , 

THEORETICAL  CONSIDERATIONS 

We  consider  the  simple  redox  reaction 
Ox  +  e  ~  s*  Rd 
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at  the  microdisk  electrode  subjected  to  an  AC  flux 

D  3c/3r-  -Q  sin  oil  say  for  0  <  r  <  a ,  r  -  0,  t  >  0  (j) 

with  zero  flux  over  the  surrounding  insulator  surface 

3c^*0  r>  a.  2  m  0,  I  >  0  (2) 

It  should  be  noted  that,  as  in  the  application  of  small  amplitude  galvanostatic  or 
potenuostauc  pulses  [1.2],  the  uniform  flux  boundary  condition  (1)  will  apply  rather 
than  any  boundary  condition  based  on  a  uniform  concentration  over  the  surface  of 
the  microdisk.  We  apply  conditions  (1)  and  (2)  to  the  solution  of  the  differential 
equation  governing  the  concentration  of  one  of  the  redox  species: 


3c  D  aV  D  3c  |  D  d2e 

3 r  dr2  +  r  dr  dz2 

(3) 

We  consider  the  f  pp1*'-*  transform  of  eqn.  (3) 

3Jc  1  3?  3J?  2-  ^  e"  A 

(4) 

where 

(5) 

s  is  the  variable  of  the  Laplace  Tranfonnatioo  and  c" 
each  of  the  redox  species.  With  the  sustitution 

is  the  bulk  concentration  of 

c-5  exp{-(ai  +  ^l)1/7r] 

(6) 

the  solution  of  eqn.  (4)  is  [l] 

*  -  ~  +  «p[  ~  («l  +  M<*)  da 

(7) 

where  the  function  g(a,q)  must  be  chosen  to  satisy  the  Laplace  transforms  of  the 

boundary  conditions  (1)  and  <2): 

3c  —  Qa  _  . 

3*  D(s2-U2) 

(81 

3c 

—  -  0  r>  a,  r  -  0 

dr 

Use  of  the  discontinuous  integrals 


and  at  ;  »  0  we  have 

- - a?*— r\(ar)Mtut) — If —  (p 

D{s^u-)J0  (a-  +  ^)'/:  ll“' 

As  in  the  analysis  of  the  chronopotentiometnc  experiment,  we  can  invert  the 
transform  before  taking  the  integral  in  eqn.  (12).  Since  we  are  interested  only  in  the 
steady  state  behavior  of  the  alternating  component  of  the  concentration,  we  confine 
attention  here  to  the  poles  i  -  ±i<a  of  the  integrand  of  the  Inversion  Theorem 

fi.\  1  /A*’*  exp(r/)dr 

/(  +  °3) 

The  residues  at  these  poles  are 

- gPl***') _  {u) 

±2iu(a4  +  ur/D2)l/*  exp(±i8/2) 
where 

tan  9  «■  u/da1  (15  j 

Thus  the  sum  of  the  residues  at  these  poles  is* 

— — — TT7J  {•  ~  »/2)l  -  i  exp[  -i(«i  -  0/2)] } 

2u(a  +ul/DI) 


u(a*  +  u2/D2)1 


sin (ut  -  9/2) 


We  apply  Cauchy’s  Theorem  to  eqn.  (13)  (i.e.  that  the  value  of  the  integrand  is  2m 
times  the  sum  of  the  residues  at  the  poles)  and  use  eqn.  (16)  in  eqn.  (12)  to  give 


r-f  ru-w&r-w* 
D  (o4  +  «V2)J)  / 


As  in  the  chronopotentiometric  case  discussed  previously,  we  evaluate  the  average 
concentration  over  the  surface  of  the  disk 

„  _1Q  fa  ,  ,  w,  ,  sin(ur  -  8/2)r  dr  da 

cAv  -  75-  /  /  Ja(ar)Jx(aa) — - '—2 — —— 

DaJ0J0  (at  +  u1/D1)W 

~2Q  (x\,,  Xi2  sin(wt  -  9/1)  da 


With  the  substitutions 
0  ~al 
1 2  -  D/u 


( 1<5| 
( 20 1 
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we  can  write  eqns.  (18)  and  (15)  as 

c  -  --q  r 

tan  6  -  1  //32  (22) 

We  can  use  eqn.  (21)  in  some  appropriate  polarization  equation  such  as  the 
Butier-Volmer  relation 


,  |  /9a  \l2sin(w/  -  9/2)  d/3 

'  n  J^r 


(21) 


Q 


(1  -  g)nr)F\ 
RT  ) 


C«+cA 


-anilFW 

RT  I} 


(23) 


where  we  have  assumed  equal  bulk  concentrations  of  the  two  redox  components. 
For  sufficiently  small  q  and  cA»  we  linearize  eqn.  (23)  in  the  usual  way  and  expand 
eqn.  (21)  into  the  in-phase  and  quadrature  components  to  give  the  real  and 
imaginary  components  of  the  AC  impedance 


-  *r  +  ART  f*f  ,/ae)12cos(g/2)d/8 

vnFi0a2  vn1F1Dl/7u'/1a1c'°  h  l  lv  1 1 J  0(1  +  fi*)1/A 

+ _ 151 _ - 

vnFiQa2  vn1F1Dx/1ux/1a1c’^^\  D  / 

**T  r«f  /  a\]»sin(»/2)dj8 

1r»JFl/>,'V'Vc«  4  / /J  /}(1+/,«)>/« 

_ 4Hr _ -  /  a2u  \ 

~  «n1F1D''Jux/1a2c~  5 1  f>  / 


(24) 


(25) 


The  functions  •«.  and  ♦,  are  Ubulated  in  Table  1  and  Fig.  1  gives  a  Cole-Cole  plot 
of  the  dimensionless  quantities 


nn1F1Dac’*‘  nFDc™ 
2  RT  Z  2  ai0 


(26) 


trrt2F2£»ac“ , 
2 RT  ( 


( azu  \ 

loia2/  ' 

i  D  I 

(22 


as  a  function  of  the  dimensionless  frequency  a2u/D.  As  expected,  the  plot  differs 
markedly  from  the  plot  for  a  large  planar  electrode  (3).  At  low  frequencies.  Z 
vanishes  as  the  transport  impedance  becomes  determined  by  the  steady  state  mass 
transfer  coefficient  to  the  microelectrode  surface.  At  sufficiently  high  frequencies 
the  results  resemble  the  familiar  plot  of  the  Warburg  impedance  and  as  u  =»  x 
Z"  =»  0  and  Z'  is  determined  by  the  charge  transfer  resistance.  This  charge  transfer 
resistance  displaces  the  Cole-Cole  plot  along  the  Z'  axis  and  R^  can  be  obtained 
from  both  the  high  frequency  and  zero  frequency  limits.  In  contrast  to  measurement 


Values  of  the  function  #.,( a:u/D)  and  ♦5(o:w/0>.  eqns.  (24M27) 


a  *u> .  D 

<Pi(a1u/D) 

♦  5(a:w/D) 

0.000400 

0.00841 

0.000070 

0.001600 

001669 

0.000276 

0  003600 

0.02483 

0.000612 

0  006400 

0.03282 

0.001075 

0.01000 

0.04068 

0,001658 

0.04000 

0.07786 

0.006213 

0.09000 

0.1116 

0.01309 

0.1600 

0.1421 

0.02179 

0.2500 

0.1695 

0.03186 

0.4900 

02156 

0.05465 

0.6400 

02341 

0.066  n 

0.8100 

0.2517 

0.07905 

1.000 

0.2665 

0.09131 

1.210 

02794 

0.1034 

1.690 

0.3002 

0.1266 

1.960 

0.3086 

0.1376 

1230 

0.3157 

0.1480 

2890 

0.3269 

0.1674 

3240 

0.3312 

0.1762 

3610 

0.3349 

0.1845 

4  840 

0.3427 

0.2064 

5.760 

0.3460 

02186 

6.760 

0.3482 

02292 

7.840 

0-3496 

0.2384 

9.000 

0.3506 

0.2464 

1024 

0.3513 

0.2534 

11.36 

0.3517 

0.2595 

12.96 

0.3520 

0.2649 

14.44 

0.3522 

02696 

16.00 

0.3524 

02739 

25.00 

0.3529 

0.2899 

36.00 

0.3532 

0.3005 

64  00 

0.3534 

0.3138 

81.00 

0.3534 

0  3182 

100.0 

0.3535 

0.3217 

1440 

0.3535 

0.3270 

1960 

0.3535 

0.3308 

256.0 

0.3535 

0.3336 

324.0 

0.3535 

0.3359 

400.0 

0.3535 

0.3376 

506.2 

0.3535 

0.3394 

625.0 

0.3535 

0.3408 

7562 

0.3535 

0.3420 

with  conventional  electrodes,  the  fequencies  are  scaled  by  the  parameter  (  D/a1 )  A' 
in  other  experiments  with  microelectrodes,  the  kinetics  of  fast  electrode  reactions 
become  measurable  by  making  (D/a)  sufficiently  large. 
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Re(Z) 

Fi tr  l-  Colt-Cote  plot  of  the  ilimninnlm  impedanrei  (24)  end  (23). 

The  equivalent  circuit  can  be  designated  by  Fig.  2  where  M  denotea  the  diffu- 
sional  impedance  of  the  microelectrode.  It  is  unlikely  that  the  uncompensated 
solution  resistance,  R„,  will  have  to  be  taken  into  account  as  the  spherical  potential 
field  in  the  solution  minimizes  Ra  (just  as  the  spherical  concenctration  field 
maximizes  km).  The  time  constant  R,  Cm  will  usually  be  short  compared  to  the 
shortest  accessible  value  of  u'1:  the  response  at  high  frequencies  will  in  fact  be 
similar  to  those  of  planar  electrodes  being  determined  by  the  product  AcrQi- 


Cdl 


i 


i 


R  cf 

Fig.  2.  Equivalent  circuit  (or  the  disk  microeiectrode 


CONCLUSION 


The  simple  approach  to  the  solution  of  the  differential  equation  governing  mass 
transfer  to  microdisk  electrodes  which  was  used  previously  to  discuss  chronopo- 
tentiometric  and  chronoamperometric  measurements  [1.2]  can  evidently  be  used  to 
analyze  all  the  usual  electrochemical  experiments  at  an  adequate  level  of  approxima¬ 
tion.  Applications  to  reaction  schemes  involving  coupled  chemical  reactions,  to 
cyclic  amperometry.  and  to  measurements  using  ring  electrodes  will  be  reported 
elsewhere  [4],  The  analysis  presented  here  opens  up  the  use  of  microelectrodes  in 
AC  impedance  experiments.  Such  measurements  are  likely  to  be  particularly  useful 
since  it  becomes  possible  to  scale  the  frequency  responses  and  charge  transfer 
resistances  by  the  parameters  ( D/a  )2  and  ( D/a )  respectively.  Therefore,  as  in  other 
applications  of  microelectrodes  to  the  measurement  of  electrochemical  kinetics,  the 
electrode  dimension  becomes  a  parameter  of  the  experimental  investigation. 
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Qi 


1 

n 

R<rr 

Z' 

Z" 
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*s 


Double  layer  capacitance.  mF 

(-D1/J 

Mass  transport  coefficient  cm  s~' 
(Z>/w)1/J 

Number  of  electrons 
Uncompensated  resistance,  0 
Charge  transfer  resistance,  Q 
Real  pan  complex  impedance 


Imaginary  pan  complex  impedance 
Tan_1(w/f>aJ) 
f*r  ,  nlJcos(«/2)  dP 


fN*?)] 


0(1  +  0*) 

2  sin(  9/2)  dP 


J o  I  v  1  'J  p(\  +  04)1 
Angular  frequency,  Hz 
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